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Abstract 

The current WRF-Fire model starts the fire from a 
given ignition point at a given time. We want to start 
the model from a given fire perimeter at a given time 
instead. However, the fuel balance and the state 
of the atmosphere depend on the history of the fire. 
The purpose of this work is to create an approximate 
artificial history of the fire based on the given fire 
perimeter and time and an approximate ignition point 
and time. Replaying the fire history then establishes 
a reasonable fuel balance and outputs heat fluxes into 
the atmospheric model, which allow the atmospheric 
circulation to develop. Then the coupled atmosphere- 
fire model takes over. In this preliminary investigation, 
the ignition times in the fire area are calculated based 
on the distance from the ignition point to the perimeter, 
assuming that the perimeter is convex or star-shaped. 
Simulation results for an ideal example show that the 
fire can continue in a natural way from the perimeter. 
Possible extensions include algorithms for more general 
perimeters and running the fire model backwards in time 
from the perimeter to create a more realistic history. The 
model used extends WRF-Fire and it is available from 

[openwf m . org] 

1. INTRODUCTION 

Fire models generally start the fire from a given 
ignition point at a given time, and sophisticated ignition 
parameterizations exist, including line ignition and 
submesh ignition procedures | |Mandel et al.| |2011| . 
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However, in practice one is often faced with the need 
to start a fire model from a fire already in progress. 
Often only perimeter data pertaining to some time 
are available, such as from the US Forest Service at 
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analyses of existing fires currently, and it wi 

even more important for forecasting of the behavior of 

fires in progress in future. 

With models that do not include two-way interaction 
with the atmosphere, continuing from an existing fire 
state is essentially straightforward. While the fuel can 
be partially burned in some areas, fuel in locations 
untouched by the fire is unchanged, and the model 
can simply progress to the new areas regardless of 
the fire history. (A model that may include long-range 
effects such as preheating in front of the fire, would 
be an exception.) In a coupled atmosphere-fire model, 
however, the situation is very different. First, simply 
igniting the whole area inside the given perimeter is not 
an option, because the large instantaneous heat release 
will cause the model to break down. More importantly, 
the state of the atmosphere evolves in interaction with 
the fire, and the buyoancy caused by the heat flux 
causes significant changes to the wind field, which in turn 
influences the future progress of the fire. 

Starting the fire from a state already developed is 
essentially a data assimilation problem, and it could be 
treated as such by a shooting method: the fire starts 
from a ignition point at a time in the past, then at the 
given simulation time, the state of the fire is compared 
with the given perimeter, and adjustments can be made 
to the ignition time and location, much as in variational 
data assimilation methods such as 4DVAR. We plan to 
study such approaches in future as a part of our effort in 
the area of data assimilation for wi ldland fires ^Beezley] 



2009 



2009 



Beezley and Mandel|r2008[ [Mandel et al.||2008 
[2010^ 



The approach adopted here is different. Given a fire 
perimeter, we create an artificial history. Then the fire 
history is replayed, which produces heat output into the 
atmosphere, and the atmospheric model spins up to 
a state that is plausible for the fire at the stage given 
by the perimeter. The artificial history is essentially 
a parameterization of the process that leads to the 
development of the fire perimeter. Parameterizations of 
various levels of sophistication can be considered, up 
to and including running a fire model backwards in time 
to find an ignition point and iterating to find a matching 
atmospheric state. In this initial study, we consider a 
very simple artificial history model, and show that it 
results in acceptable fire and atmosphere states for the 
perimeter. Approximate perimeter states obtained by 
such method could provide also a good starting point for 
data assimilation in future. 

2. THE MODEL 



Fire models range from tools based on |Rothermel|||1972| 
fire spread rate formul as, such as BehavePlus (Andrews' 
|2007j and FARSITE (Finney] |1998], suitable for oper- 



ational forecasting, to sophisticated 3-D computational 
fluid dynamics and combustion simulations suitable for 
research and reanalysis, such as FIRETEC i jLinn et alT 
|2002i and WFDS HMell et al.||2007) . BehavePlus, the 
PC-based successor of the calculator-based BEHAVE, 
determines the fire spread rate at a single point from fuel 
and environmental data; FARSITE uses the fire spread 
rate to provide a 2-D simulation on a PC; while FIRETEC 
and WFDS model combustion in 3D, which is much more 
expensive. See the survey by |Sullivan] | |2009) for a num- 
ber of other models. 

The model considered here co uples the mesoscale 
atmospheric code WRF-ARW (Ska marock et aL|[2008) 
with a fire spread module, based on the Rothermel 
model ( |Rothermel„197 2) and implemented by the level 
set method. In each time step, the fire model inputs 
the atmospheric winds and outputs surface sensible 
and latent heat fluxes into the atmosphere. Only the 
finest domain in WRF is coupled with the fire model. 
The fire model works in conjunction with WRF land 
use models, and it interpolates horizontal winds from 
the ideal logarighmic wind profile to appropriate heights 
above the surface, for each fuel. 



The model has grown out of NCAR's CAWFE code 
^Clark et al.||1996b|a| [2004|[Coen||200 5), which couples 
the Clark-Hall atmospheric model with fire spread 
implemented by tracers, and it got its start from a 
prototype code coupling the fire model in CAWFE with 
WRF in LES mode | |Patton and Coen|[2004) . The 



tracers, however, were replaced by a level set method, 
which we considered more flexible and more suitable for 
data assimilation and WRF parallel infrastructure. The 
coupled model is capable of running faster than real time 
in LES mode, with resolution of tens of meters for the 
atmosphere, and meters for the fire, with the matching 
time step of a fraction of a second, on the innermost 



modeling domain of many kilometers in size yordanov 
et al.|201 1) . Fuel data and topography can be obtained 
from government databases in the United States and 



from satellite images and GIS elsewhere i jJordanov 
erar]|20TTt . See [Mandel "eTaT] f2009 , 26TT} for futher 
details and references. The model is currently available 
from the Open Wildland Fire Modeling environment at 
l openwfm.orgl along with Utilities for data preparation, 
visualization, and diagnostics and a wiki with many user 
guides for the specific features and utilities, discussion, 
and support. A code containing a subset of the features 
is distributed with WRF as WRF-Fire. 

3. ENCODING AND REPLAYING THE FIRE HISTORY 

The state of the fire model consists of a level set function, 
$, given by its values on the nodes of the fire model 
mesh, and time of ignition t;. The level set function is 
interpolated linearly. At a given simulation time t, the fire 
area is the set of all points {x,y) where ^{t,x,y) < 0. 
The level set function and the ignition time satisfy the 
consistency condition 



^it,x,y)<0^=^Tiix,y)<t, 



(1) 



as both of these inequalities express the condition that 
the location (x, y) is burning at the time t. In every time 
step of the simulation, the level set function is advanced 
by one step of a Runge-Kutta scheme for the level set 
equation 

where R = R{t,x,y) is the fire rate of spread, which 
depends on the fuel, wind speed, and slope. The ignition 



time at nodes is tinen computed for all newly ignited 
nodes, and it satisfies the consistency condition Ij}. 

The fire history is encoded as an array of ignition times 
Ti {x,y), prescribed at all fire mesh nodes. To replay the 
fire in the period < t < Tp^, the numerical scheme for 
advancing $ and Tj is suspended, and instead the level 
set function is set to 

^{t,x,y)^Ti{x,y)~t. 

After the end of the replay is reached, the numerical 
scheme of the level set method is started from the level 
set function $ at f = Tpci-. 

For reasons of numerical accuracy and stability the 
level set function needs to have approximately uniform 
slope. For example, a very good level set function, which 
has slope equal to one, is the signed distance from a 
given closed curve r, 

$(a:;,2/) = ±dist((a;,y),r), 

where the sign is taken to be negative inside the region 
limited by r, and positive outside | |Osher and Fedkiw 
2003| . Thus, the ignition times Ti need to be given 
on the whole domain and they need to be such that Ti 
decreases with the distance from the given perimeter 
inside the fire region, and increases outside. The ignition 
times Ti outside of the given fire perimeter are perhaps 
best thought of as what the ignition times might be in 
future as the fire keeps burning. 

4. CREATING AN ARTIFICIAL FIRE HISTORY 

The purpose of this algorithm is to create the artificial 
values of the time of ignition on the fire model mesh, 
given ignition point (xign, j/ign), ignition time Tig„, fire 
perimeter r, and the time when the fire reached the 
perimeter Tpcr, assuming that the fire perimeter is 
convex, or at least star-shaped with respect to the 
ignition point. The fire perimeter is given as a set of 
points {xk,yk) in the fire model domain, k = l,...,n 
which form a closed curve consisting of line segments 
[{xk,yk),{xk+i,yk+i)] between each two successive 
points. We take ixi,yi) = {xn+i,yn+i) so that 
the starting and the ending point are identical. The 
coordinates of the point of ignition and of the points 
defining of the fire perimeter do not need to coincide with 
mesh points of the grid. 



The method consists of linear interpolation of the 
ignition time between Tign at the ignition point and Tpor 
on the perimeter, along straight lines connecting the 
ignition point with points on the perimeter. The ignition 
time is also extrapolated beyond the perimeter in the 
same manner to provide a suitable level set function, as 
discussed in the previous section. Given a mesh point 
with coordinates (x, y), the algorithm to determine the 
ignition time Ti {x,y) consists of the following steps. 

1. Find the intersection (a::b,yb) of the fire perimeter 
and the half-line starting at the ignition point and 
passing through the point {x,y) (Fig. [l]. For this 
purpose, we use the function 

F{x,y,Xb,yb) = (yfc-yign)(a:-a::ign)-(a;6-a:^ign)(y-yign), 



lies on the line 
and it is positive 



which is zero if point (a;b,2/b) 
defined by {x,y) and {xign,yign) 
in one half-plane and negative in the other. We 
then find segment [{xk,yk),{xk+i,yk+i)] such that 
F{x,y,Xk,yk)F{x,y,Xk+i,yk+i) < 0, that is, the 
points {xk,yk) and {xk+i,yk+i) lie on opposite sides 
of the line passing through {x,y) and ixign,yign)- 
Since the line intersects the fire perimeter at two 
points, one on each side of the ignition point, we 
choose correct segment as follows: 

• If {x,y) is inside r, that is, closer to the ignition 
point than to the intersection, then the desired 
segment is the one that lies on the same side 
from the ignition point as the point {x,y); 

• If {x,y) is outside of r, then the needed 
segment lies on the same side from the ignition 
point as {x,y). 

Calculate the time of ignition of the mesh point, 
based on the ratio of the distances of the mesh point 
and the perimeter point to the ignition point. 



Ti{x,y) =Tign- 



\\{x,y)~ (a:ign,yign)H 
|(a;b,yb) - (a;ign,2/ign)|| 



V-'pcr -'ignj 
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Fig. |3] we show 3D renderings of tine simulation. Tine 
streamlines near the surface show the updraft created 



Figure 1: Construction of intersection of the fire 
perimeter and the half line originating from the Ignition 
point and passing through a given mesh point. 



5. COMPUTATIONAL RESULTS 

We have tested this algorithm on an ideal example 
to measure the difference in the atmospheric winds 
between a simulation propagated naturally from a point 
and another one advanced artificially. In this example, 
the topography was flat except for a small hill roughly 
500 m in diameter and 100 m high in the center of a 
domain of size 2.4km x 2.4km. The atmospheric and 
fire grid resolutions used were 60 m and 6 m respectively, 
with a 0.25 s time step. The background winds were 
approximately 9.5 m/s traveling southwest at the lowest 
atmospheric layer 30 m above the surface. The first 
simulation was ignited from a point in the northeast 
corner of the domain 2 seconds from the start, and 
the fire perimeter was recorded after 40 minutes. This 
perimeter and ignition location were used to generate 
an artificial history for the first 40 minutes, which was 
replayed in the second simulation. Therefore, the 
fire perimeters in both simulations are identical at 40 
minutes. Both simulations were then allowed to advance 
another 28 minutes, using the standard coupled model. 
The outputs were then collected for analysis. 

Any differences in the simulations after this time are 
a result of the error of the artificial fire propagation. In 



as a result of the heat output from the fire. In Fig.|3a| the 
fire is affecting the atmosphere despite being propagated 
artificially. A semi-transparent volume rendering of 
QVAPOR was added to simulate the smoke release. 
In Fig. |2] the differences in the wind between the two 
simulations at 68 minutes and the fire perimeter are 
shown. Fig.|2a] shows the difference of the wind from the 
direct fire propagation minus the wind from the artificial 
propagation. Fig. [2bl shows the relative error in the wind 



speed defined as the norm of the difference from Fig. 2a 



divided by the wind speed from the direct simulation. 
This shows that the maximum error at the end of this 
68 minute simulation is less than 2.5%. In this case, 
the Froude number is about Fc = 0.79, showing that 
the heating from the fire may significantly affect ambient 
wind, therefore small differences caused by using the 
artificial history have an effect. The effect is concentrated 
downwind from the fire, as it could be expected. 

6. CONCLUSION 

We have presented a parameterization of the fire up 
to a point in time when a fire perimeter is given. 
The parameterization allows for the changes in the 
atmospheric circulation to develop, caused by the heat 
flux from the fire. This provides appropriate starting 
conditions for the computation to continue with the full 
coupled fire-atmosphere simulation. We have shown 
on an ideal example that the differences in the state of 
the atmosphere between a complete fire simulation and 
when the parameterization is used are not significant. 
In the studied case, the coupling between the fire and 
the atmosphere was strong. It would be interesting to 
observe how the the differences change if the problem 
moves from the wind-driven {Fc > 1) to the plume 
driven regime {Fc < 1) regime. This will be studied 
elsewhere. We plan to study also algorithms for more 
general domains, not just star-shaped, and to take into 
account different rates of fire propagation due to fuel 
nonhomogeneity. 
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